Markovian diffusion approach to the relativistic ideal gas system 
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It is shown that the usual identification of certain averages with the relativistic thermodynamical 
functions is not possible in the moving reference frame description. The Brownian motion approach 
is used as the Markovian kinetic description of the relativistic ideal gas system. 
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The number of phenomena, the critical slowing down 
[1], noise-induced transitions [2], stochastic resonance [3], 
ratchet dynamics [4], and resonant activation [5], have 
been discovered by studying the effect of the noise in 
nonlinear systems [6]. The mechanical model for the 
most of applications is the (Newtonian) Brownian parti- 
cle in a certain external potential. In spite of the signifi- 
cant progress with understanding the nonlinear stochas- 
tic behavior the rigorous results, going beyond the single- 
event stationary description, are not available for multi- 
variable systems and thus usually the overdamped ki- 
netic is considered. The stochastic description is then 
provided by the first order Langevin equation or, in the 
case of (thermal) Gaussian white noise (GWN), by the 
corresponding (Smoluchowski-) Fokker-Planck equation 
[7] in the configuration space. Note that the pioneering 
works [8] on the free Brownian motion also deal with the 
diffusion, but in the momentum space, 



Pt = -jv + €u 



(1) 



where the kinematic friction is proportional to the veloc- 
ity v = p/m and (&£ s ) = 2DS(t - s). The Eq. (1) leads 
to the correct stationary state of thermal equilibrium 



W(p) OC e -(7/-D)p 2 /2m = e - 



(3e 



(2) 



of an ideal gas system and provides the fluctuation- 
dissipation relation j/D = [3 (or D — 'yksT). The first 
order Langevin equation in the momentum space may be 
thus considered as a diffusion Markovian model of ther- 
mal equilibration for spatially homogeneous free parti- 
cles systems. The main idea of the present paper is to 
apply the same approach to the ideal gas of relativistic 
particles. The relation between momentum and velocity 
v = de/dp = pj a/p 2 + to 2 , where e = \J p 2 + m 2 is the 
energy of the particle, is then nonlinear, as well as the 
kinetic equation 



Pt 



7P/Vp 2 + m 2 + £ t . 



(3) 



The Gibbs-Boltzmann ID normalized stationary distri- 
bution easily follows from Fokker-Planck theory [7] 



W{p) 



2mK 1 (m'y/D) 



-(l/D)y/p 2 +m 2 



(4) 



and the Einstein relation j/D = (3 remains unchanged. 
K\(m(3) is a modified Bessel function. 



The use of the diffusion description in thermodynam- 
ics is related to the Markovian approximation of the ki- 
netic theory. The general Kramers-Fokker-Planck ki- 
netic equation in the phase space, see e.g. Ref. [9], sim- 
plifies to the form 



d t P(p, t) = 7 V o (pP/e) + DV 2 P 



(5) 



for a spatially homogeneous 3D system of a free rela- 
tivistic particles with a (reduced) potential interaction 
with the thermal bath. The potential interaction corre- 
sponds to the state-independent diffusion at the level of 
equation (5). V = d/dp and e = \J p 2 + m 2 . The 3D 
generalization of the Langevin description (3) is obvious, 
however in contrast to the Newtonian case it results with 
the set of three coupled equations with additive indepen- 
dent noises. The second remark concerns the specific re- 
quirements of the special relativity theory. Note that the 
Gibbs-Boltzmann equilibrium distribution is established 
in the distinguished resting reference frame of the thermal 
bath (and of the considered system as a whole), in which 
the system and the bath remain in contact and in equilib- 
rium. The parameter [3 — 1/ksT is given by the temper- 
ature of the thermal bath [10]. The any corresponding 
probability density distributions in a momentum space, 
W(p) and W'(p") in a resting and moving (along x-axis) 
with a constant velocity V reference frame, respectively, 
are related by [11] 

eW(p) = e'W'(0), (6) 

where e and p should be expressed as functions of e' and 
Jf, according to the Lorentz transformation 



(e' + Vp' x )/Vl - V 2 



(7) 



Px = (p' x + Ve')/y/l- V 2 , p y = p' y , p z = p' z . (8) 

Because the time increments, counted from the initial 
preparation of the system, are given by t — t' / VI — V 2 
and t' , respectively, the Eq. (6) leads to following relation 



1 + Vp'Je 1 

Vi - v 2 1 Vi- v 2 ' 



(EL 

VVT 



r 



Py, Pz 



(9) 

between an arbitrary P(p, t) and the corresponding 
P'(f,t'). Because the stochastic evolution is Marko- 
vian the process is completely determined in both refer- 
ence frames if the time dependent single-event probabil- 
ity density P(p, t) is known. The resulting from Eqs. (5) 
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and (9) Fokker-Planck equation in a moving reference 
frame is however very complicated in 3D, so the next 
consideration are restricted to one dimensional case. 

The (ID) equilibrium distribution in a moving refer- 
ence frame results directly from Eqs. (4), (6), (7), and 
(8), 



W'(p') 



l + Vp'/yp' + m 2 
2mVl-V 2 K 1 (ml3) 



exp 



vV 2 + m 2 + vy 



(10) 

and thus it is not involved with particular kinetic model. 
From the viewpoint of the probability theory the differ- 
ence between distribution (10) and the Boltzmann dis- 
tribution (4) is essential, i.e., it cannot be regarded by 
the change of the location parameters [12], in contrast 
to the Newtonian case. The number of averages may by 
computed from the generic formula [13] 



/dx- 
-oo 



*-bx 



\Jx 2 + m 2 



2K (my / a 2 - b 2 ), (11) 



where a > \b\, by the successive differentiation with re 
spect to a and/or b. 

-mK[(m/3) 



(£') = 



(p) 



^\-V 2 K Y (m^ 
mVK[(mp) 



(12) 



(13) 



<e' 2 > 



(A 



(pV 



/3 2 (i-y 2 )Ki 



[m 2 p 2 K'{ - m(3V 2 K[ + V 2 Kx], 

(14) 



f3 2 (l- V 2 )K 1 

_ V 

~ p 2 (l-V 2 )K 1 



[m 2 p 2 V 2 K'^-ml3K' 1 +K 1 }, (15) 



[m 2 f3 2 K '{ - mPK[ + K x \ . (16) 



The omitted arguments of (14-16) are equal m(3. The 
averages in the resting reference frame are obtained by 
putting V = in Eqs. (12-16). 

The stationary averages of twovectors 

((e,p)\ = (-mK 'jKu ) and ((e', p')) 

{(e) /Vl - V 2 , -V(e}/Vl - V 2 ), in the resting and 
moving reference frame satisfy the Eqs. (7, 8). We have 
also (p') = — V(e'). Thus the equilibrium mean energy 
and mean momentum of the Brownian particle may be 
considered as the components of appropriate twovector 
of a free particle with renormalized resting mass 

m(T) = -mK' x (m{5)lK x (m{i) [= (e)\. (17) 

The Eqs. (14) and (15) show that the invariant relation 
(e' 2 — p' 2 ) = m 2 holds identically, because the equation 
defining the Bessel function Ki(m/3), 

(mpfK'l + (mf3)K[ - [1 + (m/3) 2 ]^ = 0, (18) 



is obtained as the condition. 

The single collision with the certain "wall" of ID box 
[of a volume (length) V] in the resting reference frame 
is associated with energy-momentum transfer (e, p) — 
(e, —p) = (0, 2p) and the period between successive 
events is 2V/v = 2eV/p. Thus the pressure is 



V = (p 2 /(V^p 2 +m 2 )) = 1//3V = k B T/V. (19) 

The corresponding momentum transfer in the moving ref- 
erence frame follows from the Lorentz transformation of 
(0, 2p) and it is equal 2pj\]\ — V 2 . Simultaneously the 
period changes according to the Lorentz dilatation, so 
finally [14] 



V = - V 2 ) = fc B T/(vVl - V 2 )- (20) 

The (thermodynamical) entropy is defined in the resting 
reference frame by 

S = -k B (\n[hW{p)/V]) 

= k B {/3(e) + ln[2mV^i(m/3)//i]}, (21) 

where the Planck constant is intr oduced for dimensional 
purposes. Using U' = (e') = U/ \J\ — V 2 and 



dU' = T'dS' - V'dV' = (TdS - VdV) / \f\ - V 2 , (22) 

in view of Eq. (20) one obtains the relation T'dS' — 
TdS /VI - V 2 . Assuming T' = a(V)T and using 
Eq. (21) the relations [14] 



T' = T, and S' = S/yJl - V 2 , 



(23) 



supporting the Landsberg choice of the invariant temper- 
ature [15] are obtained. 

The identification of averages (19) and (21) with 
thermodynamical functions is not possible in a mov- 
ing reference frame. Neither the quantity defined by 
-k B (1n[hW'(p')/V']) is equal to the entropy (23) 5', nor 
the quantity (p' 2 / (V\/p' 2 + m 2 )) is equal to the pressure 
(20) V, which follows from explicit expressions. For the 
same reasons the relation between the energy fluctuation 
and the heat capacity 



Var(e) = -d(e)/d/3, 
do not hold in the moving reference frame, 



(24) 



Var(e') = 
and 



(1 - V 2 )^ 
8(e') _ m 2 K'( - (mK' x ) 2 lK x 



d/3 VI - V 2 K X 

The r.h.s. of Eqs. (25, 26) are different, if V ^ 0. 



(25) 
(26) 
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The Fokker-Planck equation of the transformed (ID) 
kinetic, Eqs. (5) and (9), reads 



c) 



"7 



d («' + V) 



P' 



+D 



dp' (1 + Vv') 2 
8 (1-V 2 ) 1/4 d (l-V 2 ) 1/4 



.dp' l + Vv> dp' l + Vv' . 
and thus corresponds to the following multiplicative 



dpf_ _ 

~d£ ~ ~ 7 ' 



v' + V (1-V 2 ) lf \ 

H ~l i T 7" / St' 



(27) 
P' 

(28) 



(i + w) 2 1 i + 

Langevin equation (written in the Stratonovich interpre- 
tation [7]), where v' = p'/^p' 2 + m 2 . The Eq. (28) 
may be obtained directly by appropriate transformation 
of stochastic processes entering to Eq. (3). We have 
d Pt = [1 + (p'/e')]dp'/Vl - V 2 and dt = dt'/Vl-V 2 , 



so 



dpt/dt = (1 + Vv')dp'/dt'. 



(29) 



Dividing the sides of (8) by the respective sides of (7) one 
has 



(v' + V)/(l + W), 



(30) 



which is the relativistic rule of collecting the vcloci- 
ti es. Be cause (£ t £ } = 2D8{t) = 2DS(t'/Vl - V 2 ) = 
VI - V 2 [2DS(t')} the GWN transforms according 



6 = a-^ 2 ) 1/4 e 



(31) 



where £t< is again a white Gaussian noise with the same 
parameter D. Collecting (29-31) we obtain the Langevin 
equation (28) in Stratonovich interpretation (because the 
ordinary calculus has been consequently used [7]). The 
Eq. (28) is a rather complicated nonlinear stochastic 
equation with (nonlincarly coupled) multiplicative (i.e., 
the state dependent) noise. For the special case of the 
ultrarelativistic kinetic (m = 0) it simplifies to the quasi- 
linear form 



dp' 



(i-y 2 ) 1/4 



dt' ~ s\gn(p') + V ' 1 + Vsign(p')^'' 



(32) 



The certain characteristics of the kinetic (28) are in 
principle known from the general Pontryagin-type equa- 
tions [16] in terms of two quadratures of the stationary 



probability density distribution and of the coefficients 
of the Langevin equation. We consider the relaxation 
time T' [17] of the stationary autocorrelation function 
(in the moving reference frame) T' = J Q dt' \(p' t ,p' ) — 

(p 1 ) 2 ] I [(p' 2 ) — (p'} 2 ] , which is a quantity known from the 
Jung and Risken formula [18]. Unfortunately, the inte- 
grals cannot be explicitly carried out for the general case 
(28), so we restrict the consideration to the ultrarelativis- 
tic kinetic (32) only. Then the calculation is elementary 
and results in 



T' = (7)8)- 1 (l - V 2 y 1/2 (5 + 8V 2 )/(2 + V 2 ). (33) 
Thus, from the viewpoint of a moving observer the decay 
of the stationary correlations proceeds slower, comparing 
to T = 5/(27/3) in the resting reference frame. Among 
the usual Lorentz dilatation some additional correction 
appears, resulting from the mentioned difference of the 
meaning of the averages in both reference frames. 



The general conclusion of the paper is the following. 
The invariance principle, in the context of the thermo- 
dynamical description of a free particles system in the 
thermal equilibrium, requires W'(p') in a moving ref- 
erence frame to be related with the Gibbs-Boltzmann 
canonical partition function W(p) (in the resting refer- 
ence frame) via Eq. (6). The resulting probability den- 
sity distribution (10) is, as the object of the probabil- 
ity theory, essentially different from the distribution (4) . 
The certain relations of the type (19), (21), (24) between 
quantities of the mechanical (mean energy, pressure) and 
the nonmechanical (entropy, temperature) origin are es- 
tablished only in a distinguished resting reference frame, 
of the considered system as a whole and of the thermal 
bath, in which they remain (in contact and) in thermal 
equilibrium. The usual identification of thermodynami- 
cal functions with certain averages is not possible in the 
moving reference frame description. The thermodynam- 
ical functions in the moving reference frame results from 
formulas (12), (20) and (23). Within the Markovian dif- 
fusion theory the kinetic is provided by Eqs. (5) and (9), 
or for the ID case explicitly by Eq. (27) or (28). The 
correlation time T' for ID ultrarelativistic case is given 
by Eq. (33). 
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